Probing Cooper pairs with Franson interferometry 



(N 

o 

(N 
C/3 



c3 



o3 



C 

o 

(N 
> 

OS 
SO 
m 

o 

(N 



X 



Vittorio Giovannetti 1 and Kazuya Yuasa 2 

1 NEST-CNR-INFM & Scuola Normale Superiore, Piazza dei Cavalieri 7, 1-56126 Pisa, Italy 
2 Department of Physics, Waseda University, Tokyo 169-8555, Japan 
(September 12, 2012) 

A setup based on the Franson optical interferometer is analyzed, which allows us to detect the 
coherence properties of Cooper pairs emerging via tunneling from a superconductor in contact 
with two one-dimensional channels. By tuning the system parameters we show that both the 
internal coherence of the emitted Cooper pairs, which is proportional to Pippard's length, and the 
de Broglie wavelength of their center-of-mass motion can be measured via current-current correlation 
measurements. 

PACS numbers: 73.23.-b, 42. 25. Hz, 73.43.-f. 74.25. F- 



I. INTRODUCTION 

Nonlocal interferometry is a highly successful tool 
for characterizing the nonclassical properties of quan- 
tum correlated systems*^ In particular, the Franson 
interferometer— suits for studying the temporal and spa- 
tial coherence of pairs of correlated particles emitted by 
a quantum source. It consists in recording time-resolved 
coincidence events at selected output ports of a couple of 
Mach-Zehnder interferometers (MZIs), each fed by one 
of the particles of each pair, whose delays are properly 
set in order to discard the first-order coherence contri- 
butions (i.e., contributions which might arise from the 
interference between different single-particle paths). In 
other words, the Franson setting is specifically designed 
to detect the global properties of the pairs (say, their de 
Broglie wavelength^—), treating them as unique objects 
which, independent of the spatial/temporal distribution 
of their constituents, can be delocalized in space and/or 
time. 

Numerous realizations of the Franson scheme have 
been implemented in optics probing either two-photon 
states produced via a parametric down-conversion 
process^— or entangled photon-hole states^ produced 
via two-photon absorbing processes. The same tech- 
nique has also been experimentally tested with pairs 
of surface plasmons excited with a two-photon entan- 
gled source^ and it has been proposed for character- 
izing coherent transport in two-dimensional mesoscopic 
conductors Interferometric setups required for such 
proposed experiments for electrons in solids have already 
been realized (e.g., see Refs. HH-H^), and several corre- 
lation experiments have been carried out^^— also with 
superconducting sources of electrons ! 32 ' 33 

Here we discuss the possibility of employing a Fran- 
son interferometer to study Cooper pairs emitted from 
a superconductor attached to a multi-port semiconduc- 
tor device. Specifically the aim of our work is to show 
how the interference fringes measured at the outputs of 
the device can be used to characterize the two-particle 
wave function of the emitted Cooper pair. Previous 
works have proposed to detect the entanglement present 



in the spin degrees of freedom of a Cooper pair via corre- 
lation measurementsi 29 i 34 ~ — The Franson interferometer 
we consider here, on the other hand, is suited for probing 
the coherence properties of the spatial degrees of freedom 
of the emitted Cooper pair. By looking at the interfer- 
ence fringes in the coincidence counts at the output ports, 
varying the relevant parameters of the Franson interfer- 
ometer, we can measure the pair correlation length of the 
emitted Cooper pair, which is proportional to Pippard's 
length characterizing the extension of the Cooper pair 
in the superconducting source, as well as the momentum 
of the center-of-mass motion of the emitted Cooper pair, 
i.e., the de Broglie wavelength of the emitted Cooper pair 
as a single object. 

This paper is organized as follows. In Sec. Hi] we re- 
call some basic facts on the Franson interferometer en- 
lightening how it allows one to characterize two-particle 
quantum coherence in the optical case for several param- 
eter configurations. Then in Sec. IIIII we discuss how it 
could be applied for the study of the Cooper-pair wave 
function. The paper ends with Sec. IIVI with conclusions 
and remarks. 



II. THE FRANSON INTERFEROMETER 

The aim of the Franson interferometric scheme is to 
reveal the coherence properties of the field emitted by 
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FIG. 1: Schematics of the Franson interferometer (optical 
implementation). The EPR-like state ([2} is injected by the 
source S along the two optical paths L and R, while the coin- 
cidence counts are recorded at the detectors, e.g., (L, — ) and 
(R, — ), to reveal the two-photon coherence of the emitted bi- 
photon state. 
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a source S which produces particles in correlated pairs, 
by recording coincidence events that take place at dis- 
tant detectors.— >^ An example of the setup is sketched in 
Fig. [TJ where S is a nonclassical source which emits two 
entangled EPR-like photons that propagate towards the 
two MZIs placed along the left and right arms of the fig- 
ure. Specifically one photon is injected along the optical 
path "L" with longitudinal wave vector fco + k and po- 
larization a, while the other is injected in the path "i?" 
with wave vector k — k and the opposite polarization a. 
Indicating the electromagnetic vacuum state with |0), 
the input of the interferometer can be expressed as 



|*) ~ |0) + |* (2) ) + --- 
with l^ 2 - 1 ) being the two-photon state 



(1) 



(2) 



where, assigning a positive momentum for a photon out- 
going from the source S, frj^/ and b^ R ^ are the 
(bosonic) annihilation and creation operators, respec- 
tively, associated with the spatial direction L/R of the 
setup. According to © the two photons of | \I> C 2 ) ^ individ- 
ually propagate with average wave vector k$, yielding a 
joint momentum p pa i r := 2hko of the pair (the de Broglic 
momentum of the system^—). Vice versa their relative 
motion is characterized by a two-particle wave function 
^fc, centered at k = with width Afc, and normalized to 
half the mode spacing 6k of the optical path, times the 
amplitude probability associated with the emitting event. 
In the space representation, this corresponds to having 
two correlated single-photon pulses with their mean po- 
sitions delocalized along the line, and with a relative dis- 
tance which instead is coherently spreading according to 
the amplitude distribution 



dk 



'2it 



(3) 



centered at x = and spread over an interval Ax ~ 1 / Ak 
(the first-order coherence length of the emitted field^— ). 

The state in (jTJ) can be seen as the state produced 
via two-photon emission by a three-level atom^ in which 
the decay of a high-energy level lvq of long lifetime yields a 
couple of correlated pulses having opposite polarizations 
and complementary energies and momenta. Alterna- 
tively |\&) can be identified with the two-mode squeezed 
state produced via a parametric down-conversion process 
through a weak but stable laser of frequency u>o = 2cfco, 
which is pumped on a short x nonlinear crystal under 
proper phase-matching conditions. In this case, in fact, 
defining the two-mode squeezing operator of the process 



fi := exp 



(4) 



with ^ k being a linear function of the pump intensity, 
the noncquilibrium steady state (NESS) of the radiation 



emerging from the crystal and propagating along the op- 
tical paths L and R reads fi|0), which reduces to ^ 
when neglecting the higher-order terms in the squeezing 
parameter vf^. For such an implementation the dereal- 
ization of the two-photon mean position corresponds to 
the coherence length of the pump (up to a few meters), 
while the typical values are Aq = 2-7r/fco — 700 nm and 



Ax ~ 500 fj,mf 



-15.44.45 



respectively. 



A. Counting coincidences 

After propagating along its corresponding interfero- 
metric arm, each photon is sent through its corresponding 
unbalanced MZI and is detected (without discriminating 
its polarization) at one of the associated output ports, 
where coincidences are recorded, say between detectors 
(L, — ) and (R, — ) of Fig. [1] (in the following we focus on 
the coincidences between these ports without loss of gen- 
erality) . The statistics of such events can be expressed in 
terms of the infinitesimal joint probability of one photon 
being revealed by detector (L, —) at time ti, = t within 
dt and the other by (R, — ) at time tn = t + r within dr, 



dpW = r/ 2 I^{x L ,t;x Rl t + T)dt 



dr. 



(5) 



where rj is the effective quantum efficiency of the detec- 
tors and 1^ (xl, t; xr, t+r) = I^ 2 \t) is the fourth- order 
correlation function, 



= — |2$(cr) cos<j>-e iip <&(cr + 6L) 

1 07T 

-e-^^cr -6L)\ 2 . (6) 

Here (■ • ■) stands for the expectation value over the input 
state |\I/( 2 )) [or |^) in (fTJ), the result being identical at 

the lowest order in the squeezing expansion], while ipcr'^ 
is the field operator at the (£ = L/R, ±) detector. Apart 
from an irrelevant contribution from the secondary input 
port of the associated MZI (see Fig. [IJ, the latter is 

V^ ±} := \[^\x h U)±e i ^{x l + 6L t ,U% (7) 

where xe and xg + 6Li are, respectively, the distances 
of the detectors from S measured along the short and 
long paths of the MZI, ipt is an extra relative phase the 
photon might accumulate while propagating through the 
MZI, and given the frequency of the mode Uk = c|fc|, 



dk 



dt) i(kx-u) k t) 



'2tt 



(8) 



is the free field along the optical path £, which connects 
S to the corresponding MZI. Finally Eq. ^ has been 
computed by setting xl = xr, 6Ll = 6Lr = 6L (we 



3 



(a) 



Q-> L 



(b) 



D 



FIG. 2: (Color online) Four interfering processes: (a) long- 
long (LL) (solid lines) and short-short (SS) (dashed lines), 
(b) long-short (LS) (solid lines) and short-long (SL) (dashed 
lines). 



assume this setting in the following), and introducing 
the parameters 



:= k 6L + (tp L + <Pr)/2, 
<P ■= (<Pl - <Pr)/2. 



(9) 



With these definitions the first term within the square 
modulus of © is identified with the interference between 
(i) a pair of photons both going the long paths of the 
MZIs and (ii) both going the short paths [LL + SS inter- 
ference; sec Fig. [3(a)]. The other two terms of © instead 
represent the processes in which one of the two photons 
goes the long path while the other goes the short path 
[LS and SL; see Fig. [3(b)]. 

In a similar way we can also evaluate the statistics of 
the single-photon detection events at the two detectors 
(L, — ) and (-R, — ) by means of the quantity 



dP^% = Vl (1) (x L/R ,t) dt, 



(10) 



which measures the infinitesimal probability of getting a 
detection at (L/R,—) in the time interval [t, t + dt] with 

i^\x L/R ,t) = j2(^ L/R '- ) Hi L/R '- ) ) 

<7 

= 4^ J dx M X ) - e^L/n+koSL)^^ + ^|2_ ( U ) 

To clarify to which extent the interference recorded in 
([6]) provides a distinctive signature of the coherent behav- 
ior of the two-photon "compound" as a unique object, we 
need to look more closely at the detection process. In par- 
ticular we need to take account of the fact that in general 
the coincidences are registered over a finite time window 
At, whose extent might be larger than the typical time- 
scale of the system (a fact that might blur the two-photon 
features). Inserting ([6]) into ([5]) and integrating the re- 
sulting expression for t € [0,T] and for t € [0, Ar], we 
compute the average number of coincidence events 



N (2) 



T r A 



dP< 2 > 



(12) 



where we assume that T is the largest time scale of the 
system (still however smaller than the coherence time 
of the pump that generates the bi-photon state). Wc 
compare this with the average number of single-detection 
events at the two detectors measured over the time in- 
terval T, i.e., 



N 



(i) 

L/R 



dP 



(i) 



L/R- 



(13) 



That is, the relevant quantity for our discussion is the 
coincidence ratio between the average number of coinci- 
dences vs the average number of couples detected by the 
two detectors in [0, T], 



K 



7V(2) 



N^N R iy 



(14) 



Under these conditions we now recognize the existence 
of several operative regimes characterized by the differ- 
ent ratios among the imbalance length SL, the first-order 
coherence length Ax of the two-photon state, and the co- 
incidence window length cAr. To analyze them we first 
fix (pn — ipL. The average number of single-detection 
events then reads 



N 



(i) 

L/R 



rjT 

47T 



dx]<f>(x) - e^Qfe + 8L)\< 



(15) 



B. Fine resolution in time Ar <C Ax/c 

Let us start by considering the ideal scenario where 
At is smaller than the first-order coherence length, i.e., 
At <C Ax/c4£ We then get a simplified expression for 



?] 2 TAt 
16tt 



|2$(O)cos0- &(8L) - <£(-SL)\ 2 , (16) 



where again the first term within the square modulus 
corresponds to the LL + SS interference while the other 
two to the LS and SL interference effects. 



1. Low imbalance regime SL <C Ai 

Suppose now that the MZIs are operated in the low im- 
balance regime SL <C Ax. This implies $(0) ~ $(±<5i), 
so that Eqs. (|16|) and ([T5"j) arc further reduced to 



AT' 



2 ^^^|4>(O)| 2 (l-cos0) 2 , (17) 
4^f^(0)(l-cos</>), 



(18) 



respectively, where we have introduced the envelope func- 
tion 



(19) 



X(8L) := / dx$*(x)$(x + 5L). 
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We notice that TV^ 2 ) exhibits oscillations in <j) of visibility 
1 and period 2w, which however can be accounted for by 
the independent single-particle interference effects taking 
places at the two MZIs. Indeed the coincidence ratio for 
the present case is given by 



K 



7AT 



(20) 



which is constant, with 7 = 7r|<f>(0)/<Y(0)| 2 being a pure 
number that depends only on the input state ([2]). 



2. High imbalance regime SL S> Ax 

The situation changes dramatically if instead the MZIs 
are operated in the high imbalance regime SL 3> Ax. 
In this case in fact, one has |$(0)| » |$(±&L)| so that 
the LS and SL contributions to can be neglected, 

while the LL + SS interference dominates the coincidence 
counts (fl"6|) , yielding 



N (2) 



ifTAr 

47T 



$(0)| 2 cos 2 ^, 



(21) 



which still exhibits oscillations in <f> of visibility 1. How- 
ever the periodicity is now n, which is half the one it 
had in the low imbalance case, and this cannot be ac- 
counted for by the single-photon events. Indeed, since 
$*(x)$>(x+6L) ~ 0, the average number of single-photon 
detections (fi~5|) is constant equal to 



N 



(i) 



r)T 



X(0), 



L ' R - 2tt 
and the ratio oscillates with period ir, 



K 



7 Ar 2 

"^COS (j 



(22) 



(23) 



Remembering © one notices that the oscillations in (|2"Tj) 
are related to the de Broglie momentum p pa i r = 2Sfco 
of the photon pair. On the other hand the fringes in 
(|17|) are related to hko (i.e., to the average momentum 
of each single component of the pair). The difference in 
the periodicity of (|21[) and (|17j) hence reflects the fact 
that in the former case the oscillations arise when the 
two-photon compound interferes with itself as a single 
object propagating with momentum 2fifco, while in the 
latter case they arise from the interference of a single 
component of the pair.— 



C. Low resolution in time At 3> Ax/c, SL/c 

Consider now the case where At is the largest time 
scale of the system (i.e., larger than Ax/c and SL/c). 
Clearly when this happens part of the information as- 
sociated with ^ is washed away. Indeed in this case 



Eq. (jT6j) gets replaced by 

-[;f(0)(l + 2cos 2 <; 



8ttc 1 



4Re X(6L) cos (/> + Re X(26L)], (24) 



where X(6L) is the envelope function defined in (fTT)| . 
In the low imbalance regime SL <C Ax, this gives 



Aire y M 



(25) 



which as (|17jl presents only single-particle oscillations, 
with Nj*j R being expressed as in (fl~8"|) . and the coinci- 
dence ratio reads 



K = 7 = ir/cTX(0). 



(26) 



On the contrary in the high imbalance regime SL 3> Ax, 
we get 



8ttc 



X(0)(1 + 2 cos 2 



(27) 



which as in (|21[) exhibits oscillations of period n but with 
a reduced visibility, with Nj^J R being expressed as in (|22l) . 
The coincidence ratio in this case is given by 



K = 4(1 + 2cos 2 < 



(28) 



which oscillates with period n, but with a reduced visi- 
bility. 

As extensively discussed in Refs. HHH this prevents 
one from detecting the coherence of the superposition ([2]) 
(indeed there exists classically correlated input states, in 
which the coherence of l^ 2 )) is broken that yields the 
same signal). Still the 7T-periodicity of the oscillations of 
(f2"T| bears a clear signature of the de Broglie momentum 

Ppair ■ 



D. Recovering the two- particle de Broglie 
wavelength in the low imbalance regime SL -C Ax 

The analysis presented in the previous section shows 
that, independent of the relative width of Ar, the worst- 
case scenario for recovering the two-photon features of 
the input state is the low imbalance regime SL <C Ax. 
Under this circumstance in fact the oscillation of 
can be just interpreted as due to the single-particle in- 
terference effect recorded in N^J R . Enforcing the high 
imbalance condition SL 3> Ax, on the other hand, sup- 
presses the single-particle interference effect while keep- 
ing the two-photon one. 

Unfortunately for some of the applications we have in 
mind (see the next section), choosing SL Ax could 
be a challenging constraint. A possible way out is to 
add some correlated noise in the interferometer which 
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suppresses the single-particle coherence effect, while pre- 
serving the two-photon coherence. For instance this can 
be done by keeping ^ + <Pr constant while randomizing 
<P = {(fL — ¥>fl)/2- By doing so in fact N^} R loses the de- 
pendence upon <j) and reduces to the constant value given 
in (|2"2"|) . while Eq. ^ gets replaced by 



2 cos 2 < 



+ |$(cr + 5L)\ 2 

+ \<S>{ct~8L)\ 2 ). (29) 



Consequently, even if 6L <C Ax, we get 
r) 2 TAr 



K 



8tt 
7Ar 

" ~2T~ 



$(0)| 2 (1 
(1 



2 cos 2 



2 cos 2 



(30) 
(31) 



for At <?C Ax/c, while Eqs. (|2"7|) and (|2"5|) are reproduced 
for At 3> Ax/c: in both cases the coincidence ratio os- 
cillates as a function of with period tt but with half 
visibility. 



III. COOPER-PAIR EMISSION 

Having clarified the basic idea of the proposal we now 
show how it could be applied for characterizing the two- 
particle wave function of a Cooper pair emitted from 
a superconductor into external metallic contacts via a 
tunneling process. The fundamental observation here is 
that, within the validity of the self-consistent mean-field 
approximation^ the BCS state |BCS), which describes 
the electronic state in a superconductor, is a fermionic 
version of the two- mode squeezed state ([1]). The main 
difference is that now the bosonic operators b k „ en- 
tering in the squeezing operator in ((4]) get replaced by 
fermionic counterparts. Specifically we can write 



|BCS) =O|0), 



(32) 



where now 



O := exp 



J d 3 kG k (e i9 al t al kii -e- 



(33) 

with ak,a and a k being the annihilation and creation 
operators for the electrons in the superconductor char- 
acterized by a spin a and momentum k, and we have 
written the phase of the superconductor 8 explicitly. At 
the first order in the tunneling process when the electrons 
leak out of the superconductor, the state |BCS) produces 
a correlated two-particle state whose structure reminds 
us of the photonic state ([5]) . We wish to detect it via the 
Franson interferometery. 



A. Setting 

The setting we have in mind assumes that two one- 
dimensional metallic leads L and R are attached to a 



superconducting source of electrons S, described by the 
BCS state {35) with and MZIs are set up on both 
sides, as depicted in Fig. [TJ 

We describe the emission of the electrons from S to L 
and R dynamically by a Hamiltonian H = H$ + Hl + 
H R + H T 

H S =j^(au ( a r) ,(34) 

H t = ]T I "dke k b^b k e l (I = L,R), (35) 

CT =t4 

n t n> A* i « v 



(36) 



where now b k and b k J are the fermionic operators for 
the electrons propagating in lead £ = L,R before the 
MZIs. A = |A|e ie is the gap parameter of the supercon- 
ductor S, and 



h 2 k 2 
2m 



Ms 



(37) 



is the energy of an electron measured relative to the Fermi 
level [is of the superconductor S. The Hamiltonian H$ 
can be diagonalized by the Bogoliubov transformation 



-fe,4-y 




u k -v k \ j Qffc,t 



(38) 



within 



JO 



—Uf-fc (39) 



and 



(40) 



The BCS state |BCS) can then be identified with the 
vacuum state annihilated by the quasiparticle operators 
ctfc.CT- Notice that the parameters Uk and v k introduced 
here are related to G k of (|33[) via the transformation 



uk = cosG/c, 



Vk 



e ie sinG fe . 



(41) 



The electrons are transferred between the source and the 
leads by the transmission Hamiltonian Ht, which anni- 
hilates an electron in the superconductor S and creates 
an electron in a lead £ = L,R, and vice versa for the re- 
verse process. We assume that the electrons are emitted 
in the same manner to the left and right leads, namely, 
the transmission matrix elements for leads L and R are 
given, respectively, by 



rpiL) rp i( — kn—k')-nd/2 



rp(R) _ rp -i(kn-k')-nd/2 
± k,k' — 1 k,k'c , 



(42) 
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where c! is the distance between the contacts with the 
two leads L and R, which are assumed to be extended in 
the opposite directions — n and n, respectively. 

The emission of the electrons starts from the initial 
state |BCS)s ® |FS) p ® |FS) r, where the electron source 
S is in the BCS state |BCS)s while the leads L and R are 
in the normal states |FS)x,//j, filled with electrons up to 
the Fermi levels Hl/Ri an d the system eventually reaches 
a NESS in the long-time limit t — > oo. We look at the 
correlations between the detectors after the MZIs on the 
L and R arms in this NESS. 



B. Correlation functions and state 

We tune the Fermi levels of the two leads fip i p within 
the gap of the quasiparticlc spectrum of the supercon- 



ducting state, i.e., 



< eV L = eV R = eV< |A|, (43) 
eVt = li S -H (£ = L,R). (44) 



In this way, the quasiparticlc emission is blocked by the 
Fermi seas of the leads L and R, and the Cooper-pair 
emission is the only allowed process. Under this condi- 
tion, collecting the contributions up to the second order 
in the transmission Hamiltonian Hp, the only two-point 
correlation functions in the NESS {t\,t 2 — > oo) that are 
nonvanishing for large xi,x 2 &re&~— 



(4 ei \x 1 ,t 1 )4" 2 \x 2 ,t 2 )) = -{4 ei \x u t 1 )4 e2 \x 2 ,t 2 )) 



im 



eV 



eV 



T (l! 



Tr 



dE / d 3 ku k v k - 



k(E)-k k(-E),h 



k(E)k{-E) E 2 - h 2 ujl + i0+ 



e z[k(E)x 1 +k(-E)x2] e -iE(t 1 -t 2 )/h 



(45) 



r 



and their complex conjugates, where ip)p (x, t) is the field 
operator for the electrons before the MZI on arm £ (= 
L,R) in the Heiscnbcrg representation, 



^\x,t)=e im ' h ^\x)e- im ' n , 



dk 



lW Jkx 



2ir 



(46) 
(47) 



which is the counterpart of © but with b*p} a be- 
ing now the fermionic operator of (|35[) , and k(E) = 
kpy/l + E/fis with the Fermi wave vector kp = 
2n/\p = yj2mns /h 2 . Under the assumptions kpd 1 
and \A\/ns < 1, we get 



(4 l) (x 1 ,t 1 )^\x 2 ,t 2 )) 



2tt 



X 2 



V F (ti 



h)), 

(48) 



where 

§(x) = & L ' R \x) = & R > L \x) 



A f eV/\A\ 



ds 



2 < J-eV/\A\ VT^ 

$(x) = $ (L < L) (x) = & H < R \x) 



eV/\A\ 



ds 



27t£ J- eV /\A\ VI - .S 2 



(49) 

(50) 



with dimcnsionless constants A = e i6 \2ir) b / 2 (m 2 /H 4 ) 
x (T 2 kF , kFn e^ d - T 2 ^_ kFn e-^ d )/k F d and A = 
i{2KYl 2 {m 2 /h 4 )$d 2 hT kF , kFn T kF 
length 



-kpn 



, and Pippard's 



(51) 



characterizing the extension of a Cooper pair^ with 
vp = hkp/m being the Fermi velocity. In the follow- 
ing, we assume d/n^ <C 1, which is typically the case: 
see Ref. HU and the caption of Fig. [31 Then, in the zero- 
bias regime eV <C |A|, the pair wave functions behave 
as 



$(x)/A~ $(x)/A 



eV 
<\A 



sinc[(xM)(ey/|A|)] (52) 



for x/7r£ 3> 1, which as a function of x has a width 
Ax ~ n 2 £\A\/eV = nhvp/eV. For eV = |A|, on the 
other hand, they are estimated to be 

${x)/A ~ <l(x)/i ~ ijo(xM) (53) 

for the whole range of x, where Jo(x) is a Bessel function. 
This has a width Ax ~ 2.4 7r£. See Fig.[3j where $(x) is 
shown as a function of x and eV. 

Notice that (ipi e) (x , t)) = O^tM and recall that the 
correlation functions in (|45l) and their complex conju- 
gates are the only nonvanishing two-point correlation 
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FIG. 3: (Color online) The pair wave function |$(a;)| 2 [nor- 
malized by (j4eV/|A|) 2 ] as a function of x and eV, on the 
basis of the numerical integration of (|49p . The parameters 
are |A| = 150 fxeV, \ F = 2-ir/k F = 0.36 nm, (i s = 11.6 eV, 
d = 150 nm, and £ = 2.8 (im for an Al superconductor.— 

functions up to the second order in the transmission 
Hamiltonian Ht- This implies that the state of the elec- 
trons in the leads L and R emitted from the supercon- 
ductor S is expressed, up to the second order in Ht, as 
O withS 

|^)):=E/**F 2) «4«,tl0)» (54) 
ti,e 3 

where ij/j^ 1 '^ are the Fourier transforms of ^ il,£ ^(x) as 
in ([3]). This should be compared with ([5J for the optical 
case, where the polarization state of a pair of photons 
is symmetric, while in the superconductor case in (|54[) 
the spin state of a pair of electrons is antisymmetric (the 
singlet state). On the other hand, the states of their 
spatial degrees of freedom arc symmetric in both cases. 

C. Coincidences 

We can now investigate the Franson interferometry for 
the electrons emitted from a superconductor. Using ((TJ) 
the coincidence rate between detectors (L,—) and (R, — ) 
is found to be ruled by an expression which reminds us 
of © for the optical scenario, i.e., 

I (2) (t) = — \2$(v f t) coa<p - e iv, <Z>(v F T + SL) 

167T 

- e - llfi <S>(v F T-SL)\ 2 , (55) 

where now 

4> := kpdL + (<p L + <p R )/2 (56) 

instead of ©. The phases <Pl/r can be induced by the 
Aharonov-Bohm effect, by applying magnetic field per- 
pendicular to the MZIsi21"^S§r^i It is clear that the for- 
mulas for the optical case are translated to those for the 



electric case by simply substituting fco — > kp and c — > vp, 
and the pair wave function $(x) of the emitted Cooper 
pair is given by (|49|) . The single-electron detection rate 
at (L/R, — ), on the other hand, is proportional to 

I {l) {x L/R ,t) 

= A / dx |$(x) - e ^L /R +k F SL)^^ x + SL ^2 (57) 

with A=1 + |A/A| 2 , which is to be compared with (fTT| 
for the optical scenario. The correction proportional to 
\A/A\ 2 is due to the possible emission of a Cooper pair 
into the same lead, cither L or R, which was assumed to 
be absent when we considered the optical scenario. Such 
processes would be suppressed by emp loying the Cooper- 
pair splitter , 40 i 54 realized in Refs. |32| andl33l. 

It is now clear that the formulas for the coincidence 
rate listed in Sec. [II] can be applied to the present electric 
scenario by the small changes fco — > kp and c — > Vp, and 
the introduction of the factor A. Then, by employing the 
strategies outlined in Sees. IIIBHIIDl we can probe the 
coherence properties of the emitted Cooper pair. That 
is: 

(i) We perform the correlation measurements, count- 
ing the coincidences between output ports {L, — ) 
and (i?, — ) of the Franson interferometer, in the 
low imbalance regime 5L <C Ax, with ipi, = ifR. In 
this case, the average number of coincidences N^> 
oscillates function of <f> with period 27r. See 
(|17l) for fine time resolution At <§C Ax/c and ([25)1 
for low time resolution At 3> Ax/c. This oscil- 
lation with period 2-k is due to the single-electron 
interference through each MZI. 

(ii) We then move to the high imbalance regime 5L 3> 
Ax, with ifL = (fa, and observe that N^ 2 ' oscil- 
lates as a function of <fi with period ir. Sec (p?Tj) 
for fine time resolution At <C Ax/c and (|27p fol- 
low time resolution At ^> Ax/c. This oscilla- 
tion with period 7r in turn is due to the interfer- 
ence between the LL and SS processes [Fig. [2{a)] , 
and it can be regarded as the MZ interference of a 
Cooper pair as a single object, as if a Cooper pair 
as a whole goes through the long and short paths 
of a MZI and interferes, even though it is actu- 
ally split into two different branches and its con- 
stituent electrons go through different MZIs. From 
the interference fringes in 

AT(2) 

as a function of 5L, 
we can measure the de Broglie wavelength of the 
center-of-mass motion of the emitted Cooper pair, 
Apair = 27r?i/pp a i r = Ti/kp = Xp/2, i.e., the de 
Broglie wavelength of the Cooper pair as a single 
(composite) particle. 

(iii) In addition, by observing the transition between 
the low and high imbalance regimes, we can mea- 
sure the pair correlation length Ax of the emitted 
Cooper pair. Since Ax is proportional to Pippard's 
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length £ [see (|49|) . ([52]) . and ([53])]. we can also esti- 
mate Pippard's length £ characterizing the size of 
a Cooper pair in the superconducting source. 

In these ways, we can probe the coherences of both the in- 
ternal and center-of-mass degrees of freedom of the emit- 
ted Cooper pair. 

It could be, however, a bit challenging to enforce the 
high imbalance condition SL 3> Ax in practice. The 
latter would be typically of the order of Ax ~ 10 /iin, 
whereas the wavelengths of the electrons passing through 
the MZIs would be Xf ~ 0.5 nm (see Fig. [3] and its cap- 
tion). It would be difficult to take such a large path 
difference SL for each MZI. Even in such a case, we 
can employ the randomization procedure outlined in Sec. 
ED] 

(iv) We perform the correlation measurements in the 
low imbalance regime SL <C Ax, with ipL + VR kept 
constant but tp = (ipi, — <Pr)/2 randomized. Then, 
iV^ 2 ) behaves as ([50]) when the time resolution is 
fine, At <C Ax/c, while it behaves as (|27p when 
the time resolution is low, At 3> Ax/c. In both 
cases, oscillates function of <p with period 
7r (with half visibility), and we can probe the de 
Broglie wavelength of the emitted Cooper pair. The 
estimation of Pippard's length £ is problematic, but 
it is at least in principle possible by studying the 
variation of the visibility of the oscillations of 
as a function of the coincidence time window At. 
Indeed from ([27]) and ([50]) we notice that, while for 
At <C Ax/c the visibility of (or equivalently 
the visibility of the ratio 1Z) is increasing with At, 
for At ^> Ax/c it saturates to a constant value. 
Monitoring this change of behavior enables one to 
estimate Ax, i.e., Pippard's length £. 

As for the spin degrees of freedom of the Cooper pair, 
the coincidence ratio 1Z defined in ([T4"]) exceeds unity, for 
4> = 0, in the ideal situation where the time resolutions of 
the detectors are fine and the accumulation time duration 
T is short. This is bunching and is due to the singlet spin 
correlation of the Cooper pair. Therefore, by observing 
the bunching, we can detect the spin entanglement of the 
Cooper pairi 29 ' 34 " — For this purpose, however, the Fran- 
son setting is not essential, and in addition the bunching 
is lost if the time resolutions of the detectors are not very 
fine and the accumulation time duration T is relatively 
long, as considered in the present paper. 

IV. CONCLUSIONS 

In this paper we reviewed some of the features of the 
Franson interferometric scheme^ discussing how it can 



be used to discriminate the two-particle coherence effects 
from the single-particle coherence effects in different ex- 
perimental settings. Subsequently, exploiting a formal 
equivalence between the state of a Cooper pair emitted 
via tunneling from a superconductor into external metal- 
lic contacts and the two-mode squeezed state generated 
in a parametric process, we have analyzed an application 
of this setup for electronic sources. 

As already mentioned the emission of the electrons 
comprising a Cooper pair into two separate metallic leads 
has been realized recentl y 32 ' 33 The main technical chal- 
lenge for the implementation of the proposed scheme is 
posed by the realization of the beam splitter transforma- 
tions and by the ability of controlling the relative phase 
accumulated in the propagation of the electrons through 
the MZIs. A beam splitter for electrons was realized in 
Refs. and[26| in a GaAs high-mobility two-dimensional 
electron gas system, where two input and two output 
ports arc defined by applying negative bias voltages to 
gates that deplete underlying electrons, and the split- 
ting ratio of the beam splitter is adjusted by applying an 
appropriate small bias voltage to another thin gate. A 
correlation experiment, which is an important ingredient 
for the Franson interferometry, was carried out in such 
a system . 25 ' 26 However, for this type of beam splitter, 
the back reflection of electrons to the input ports is in- 
evitable, resulting in a reduction of the number of coinci- 
dence events. This problem is absent in the beam splitter 
realized for the chiral edge current of a two-dimensional 
electron gas in the quantum Hall effect regime^L r 24 ' 27 ' 28 
The quasi-one-dimcnsionality of the edge current is also 
appropriate for constructing interferometers, and MZIs 
were actually realized 2 ^— and correlation experiments 
were also performe d 27 ' 28 in such systems. Furthermore 
the possibility of interfacing quantum Hall bars with su- 
perconductors at hi gh mag netic fields has been recently 
established in Refs. [55l - [58l Using chiral edge currents 
instead of metallic leads appears hence to be a possible 
option for the implementation of our setup. 
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